Non-ohmic spin transport in n-type doped silicon 
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We demonstrate the injection and transport of spin-polarized electrons through n-type doped sili- 
con with in-plane spin- valve and perpendicular magnetic field spin precession and dephasing ( "Hanle 
effect") measurements. A voltage applied across the transport layer is used to vary the confinement 
potential caused by conduction band-bending and control the dominant transport mechanism be- 
tween drift and diffusion. By modeling transport in this device with a Monte-Carlo scheme, we 
simulate the observed spin polarization and Hanle features, showing that the average transit time 
across the short Si transport layer can be controlled over 4 orders of magnitude with applied voltage. 
As a result, this modeling allows inference of a long electron spin lifetime, despite the short transit 
length. 

PACS numbers: 85.75.-d, 72.25.Dc, 72.25.Hg, 85.30.Tv. 
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It has been a longstanding goal in semiconductor spin- 
tronics to inject, transport, manipulate, and detect spin- 
polarized carriers in silicon-based devices.— i^i^ Despite 
great success in the field over the past ten years using 
other semiconductors such as GaAs^i^i^i^iiSii^ the goal of 
achieving the same with Si has been reached only re- 
cently, using all-electrical hot-electron methods with un- 
doped single-crystal silicon transport layers. Later, spin 
injection into silicon was realized as well in an epitaxially- 
grown silicon n-i-p diode structure using circular polar- 
ization analysis of weak electroluminescence spectra for 
spin detection across a transport layer of 80-140nmjii 
Although our previous studies demonstrate electron spin 
manipulation in undoped silicon - even over a very long 
distance (350 microns )^^'^'^i^^ - it is necessary to investi- 
gate magnetic- and electric-field control of electron spin 
in doped silicon for integration of spintronics into present- 
day silicon-based microelectronic technology, where im- 
purity doping plays a critical role. 



In this report, we present spin injection, transport and 
detection in an n-type doped silicon device using our all- 
electrical methods. Unlike previous studies with undoped 
Si, the presence of ionized impurities in the depletion re- 
gions of these doped transport layers gives rise to conduc- 
tion band bending that for sufficient biasing conditions 
confines injected electrons for long dwell times. By mod- 
eling transport with drift and diffusion in the inhomoge- 
neous electric fields provided by the band bending with 
a Monte-Carlo method, we simulate both spin precession 
and spin decay, showing that the transit time distribu- 
tion of spin-polarized electrons can be controlled over a 
very wide range with an applied voltage, and can yield a 
measurement of spin lifetime. 
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FIG. 1: Schematic side view of n-type doped Si spin transport 
device and illustration of the components and contacts for 
electrical injection and detection. 



II. EXPERIMENTAL METHODS 

Fig. [T] illustrates the structure of our device. Fabrica- 
tion consists of ultra-high vacuum metal film wafer bond- 
ing to assemble a semiconductor-metal-semiconductor 
hot-electron spin detector; a silicon-on insulator (SOI) 
wafer including a 3.3/im single-crystal (100) nominally 
1-20 ri-cm phosphorus-doped n-type sihcon spin trans- 
port layer is bonded to an n-type bulk silicon collector 
wafer with a Ni8oFe2o (4nm)/ Cu (4nm) bilayer. Con- 
ventional wet-etching techniques expose the SOI device 
layer, onto which a ferromagnetic-emitter tunnel junction 
hot-electron spin injector is built. The final device struc- 
ture is Al (40nm)/Co84Fei6 (10nm)/Al2O3/Al (5nm)/Cu 
(5nm)/n-Si (3.3 /im)/Ni8oFe2o (4nm)/Cu (4nm)/n-Si 
substrate, as displayed in Fig. 1. Further details on 
fabrication of similar devices can be found in previous 
report siSii^ii^. 

An applied emitter voltage Ve on the tunnel junc- 
tion (larger than the Cu/n-Si injector Schottky barrier) 
injects hot electrons tunneling from the ferromagnetic 
Co84Fei6 cathode through the thin-film Al/Cu anode 
base and into the doped silicon transport layer conduc- 
tion band. The first collector voltage (Vci) controls 
the voltage drop across the transport layer and modifies 
the spatially nonlinear conduction band potential energy. 
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Electrons escaping the transport layer are ejected over 
a Schottky barrier at the detector side into hot-electron 
states in a buried Ni8oFe2o thin film. The final spin polar- 
ization is detected by measuring the ballistic component 
of this hot electron current (second collector current, Ic'2) 
in the n-type Si wafer below; spin-dependent scattering 
in the ferromagnetic NigoFe2o makes this current depen- 
dent on the projection of final spin angle on the Ni8oFe2o 
detector magnetization. 



III. EXPERIMENTAL RESULTS 

The spin-detection current Ic2 was first measured with 
an external magnetic field parallel to the device plane. A 
spin- valve effect, resulting from the different in-plane co- 
ercive fields of injector and detector ferromagnetic layers, 
is displayed in Fig. [2l The measurements were done with 
Ve = -1.6V applied, using different values of Vci between 
4.5V and 8V at temperature T — 152K. Because of the 
Ic2-Vci dependence, we normalize the data for compar- 
ison between different Vci values. After this normaliza- 
tion, it can be seen that the measurement is only weakly 
dependent on accelerating voltage Vci over this range. 

The in-plane magnetic field was swept between -4 kOe 
to -f 4 kOe for this measurement. Since the coercive fields 
of both ferromagnetic (FM) layers are smaller than 200 
Oe, the data obtained from the Vci = 5V measurement 
is magnified in the inset of Fig. 2 and the field sweep 
direction is specified by correspondingly colored arrows. 
When the in-plane magnetic field reaches approximately 
-f 20 Oe from the negative saturation field (below -300 
Oe), the NigoFe2o layer switches its magnetization, caus- 
ing an anti-parallel (AP) configuration in the two FM 
layers, which lowers the Ic2 current relative to a parallel 
(P) configuration, because in this case spin "up" is in- 
jected, but spin "down" is detected. If the magnetic field 
increases further, the Co84Feig layer reverses magnetiza- 
tion, resulting in a P configuration and restoration of the 
higher Ic2- This happens as well in the opposite sweep- 
ing field direction due to the symmetric but hysteretic 
coercive fields of each FM layer. The magnetocurrent 
(MC) ratio [Ic2^- Ic2^^)/Ic2^^ calculated from the 
spin- valve plot, where the superscripts refer to P and AP 
magnetization configurations in the two FM layers, is ap- 
proximately 6%. As the magnetic field reaches up to ±4 
kOe after the magnetization reversal of both FM layers, 
Ic2 monotonically rises because of domain magnetization 
saturation in the direction of the external field. 

To unambiguously confirm spin transport through 
the doped silicon layer, we have performed measure- 
ments of Ic2 in an external magnetic field perpendic- 
ular to the device plane, which allows us to exam- 
ine spin precession and dephasing (Hanle effect) dur- 
ing transport .■ia>i^>i^>i^>i^>i^>^>i^>i^>20 Depending on the 

magnitude of the applied magnetic field and the tran- 
sit time (subject to drift and diffusion through the con- 
duction band from injector to detector), the polarized 
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FIG. 2; Spin valve effect in doped Si spin-transport devices. 
Measurement was done at T = 152K and Vb = -1.6V, with 
different Vci values applied as indicated in the plot. Inset; 
Data measured with Vci = 5V plotted over a smaller field 
range. Field sweep directions are indicated by red (increasing) 
and blue (decreasing) arrows. 



electron spin (initially parallel to the injector FM layer 
magnetization) can arrive at the detector having rotated 
through precession angle d = rg^sB /Ti^ where t is the 
transit time, B is the magnetic field, g is the electron spin 
g- factor, is the Bohr magneton, and h is the reduced 
Planck constant. 

Our measurements in a perpendicular magnetic field, 
using the same experimental conditions as were applied 
in the spin- valve effect measurement (Ve = -1.6V and T 
= 152K), are shown in Fig. [3]for the same varied values 
of Vci as in Fig. [51 The measured Ic2 was normalized for 
data comparison at different accelerating voltages Vci, as 
in the spin-valve effect experiment. Again, the inset of 
Fig. [3] shows the data for Vci = 5V with magnetic field 
sweep directions indicated by correspondingly colored ar- 
rows. When a perpendicular magnetic field sweeps from 
-4 kOe (or from -1-4 kOe), /c2 exhibits a minimum be- 
fore the field reaches Oe and then it suddenly drops and 
slowly moves up between and -|-lkOe. The former min- 
ima is induced by a full spin flip due to spin precession 
(average tt rad rotation) during transport through the 
doped silicon layer, and the latter is induced by the in- 
plane magnetization switching of the two FM layers by a 
residual in-plane component of the largely perpendicular 
magnetic field, causing an antiparallel injector/detector 
magnetization configuration and reduction in signal as 
seen in previously-discussed in-plane spin- valve measure- 
ments. This argument is further upheld by changing Vci; 
minima attributed to precession appear at higher mag- 
nitude of applied perpendicular magnetic field as Vci in- 
creases due to the shorter transit time, while the FM 
switching fields clearly do not change 
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FIG. 3: Electron spin precession in doped Si spin-transport 
devices. Emitter voltage and temperature are same as in Fig. 
[21 Minima, corresponding to tt rad precession angle appear 
at higher magnetic field as Vci increases and transit time de- 
creases. Inset: data measured with Vci = 5V. Field sweep 
directions are indicated by red (increasing) and blue (decreas- 
ing). 



IV. EXPERIMENTAL ANALYSIS 

Average spin transit times r on the order of 45 - 180ps 
can be determined from the magnetic field values at tt 
rad precession minima in Fig. 3 (~ I kOe - 4 kOe) 
using^" T — h/2gfiBBTr, where h is the Planck constant. 
Correlating spin polarization from spin-valve measure- 
ment to these transit times can, in principle, be used 
to determine spin lifetime. However, these transit times 
are very short, so direct correlation as in Rcf.^^ is un- 
able to independently determine the (long) spin lifetime 
of conduction electrons in doped Si. We have previously 
measured spin lifetime of 73ns at similar temperature us- 
ing a 350 micron-thick undoped Si transport layer device; 
this lifetime increases to over 500ns at 60KJ^ 

In the undoped silicon transport layers used in previ- 
ous works jiiSiiSiii^iiiiiSS the Schottky depletion region was 
much larger than the layer thickness. Therefore, the 
conduction band was linear, resulting in a spatially con- 
stant induced electric field, and relatively "ohmic" spin 
transport where the spin transit time was inversely pro- 
portional to the injector-detector voltage drop. In this 
work, however, carrier depletion of the doped silicon due 
to Schottky contacts and the resulting space-charge from 
ionized impurities causes a nonlinear conduction band 
that can have a potential energy minimum between de- 
pletion regions unless the voltage drop is very large. Since 
injected electrons may sit in this potential well for a long 
time before escaping over the detector barrier, their spins 
will depolarize and the observed MC ratio will be sup- 



FIG. 4: Depletion-approximation conduction-band diagrams 
of the doped Si spin transport layer with injector-detector 
voltage drop of (a) 3.0V (where the transport is dominated 
by diffusion against the electric field at the detector side); 
(b) 4.46V (where the bias is enough to eliminate the neutral 
region and fully deplete the transport layer); and (b) 6.0V 
(where the potential well has been eliminated and transport 
is dominated by drift in the unipolar electric field). 



pressed. 

To significantly reduce this dwell time, an accelerating 
voltage (induced by applied voltage Vci, which adds to 
approximately 0.3V of the applied emitter voltage due 
to resistive tunnel junction electrodes^") can be used to 
alter the confining potential energy. In particular, for suf- 
ficient voltage the confining potential can be eliminated. 
It is therefore expected that the spin signal is strongly 
sensitive to applied voltage and "non-ohmic" spin trans- 
port results. 



V. MODEL 

Modeling this non-ohmic behavior is necessary. In pre- 
vious works using undoped Si transport layers where the 
electric field is constant from injector to detector, a mod- 
eling technique using the arrival-time distribution given 
by the Green's function solution to the drift-diffusion 
equation can be easily implemented. ^^i^^'^^i^^ However, 
the electric field in these doped Si devices is highly inho- 
mogeneous, making it difficult to implement the standard 
method here because the drift velocity is spatially depen- 
dent, requiring Green's function solution of a nonlinear 
partial differential equation. In general, this procedure 
is non-trivial. 

To overcome this problem and simulate spin transport 
behavior in these doped devices, we use a Monte-Carlo 
technique which translates electrons a distance v{x)At 
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FIG. 5: Monte-Carlo simulated transit-time distributions for 
injector-detector voltage drops of (a) 3.0V; (b) 4.46V; and (c) 
6.0V. Note timescale changes over 4 orders of magnitude from 
(a) to (c). 



(due to drift), and zL\/2DAt (due to diffusion) in a 
timestep At, where v{x) is tlie drift velocity at the po- 
sition X and D is the diffusion constant. (The sign on 
the latter expression is randomly chosen to simulate the 
stochastic nature of 1-dimensional diffusion.) 

The spatially-dependent electric field is calculated 
within the depletion approximation. Using a doping den- 
sity of 7.2 X 10^'* cm~^, injector Schottky barrier height 
of 0.6eV (for Cu/Si) and detector Schottky barrier height 
of 0.75eV (for NiFe/Si) results in a band diagram whose 
dependence on injector-detector voltage drop is shown 
in Fig. m This figure illustrates that the voltage drop 
across the Si transport layer can be used to alter the 
dominant transport mode: at low bias a wide neutral re- 
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FIG. 6: Results of Hanle effect simulation, Eq. ([T]). Similar to 
the experimental results in Fig. [S] minima corresponding to 
TT rad precession angle appear at higher magnetic field as Vci 
increases and transit time decreases. Magnetization switching 
is simulated by signal sign reversal between the coercive field 
values of the injector and detector ferromagnets. 



gion exists between depletion regions and electrons must 
diffuse against an electric field to escape to the detector, 
whereas for biases greater than 6V, the potential mini- 
mum is annihilated by the boundary so the internal elec- 
tric field carries electrons toward the detector everywhere 
and drift is expected to dominate. 

A realistic empirical mobility model using Eq. 10 
from Refi^ is used to evaluate v{E). The diffusion co- 
efficient at each point in space is then calculated from 
the Einstein relation D{x) = ii{x)kBT/q, where mobility 
/i(x) = v{x)/E{x) and E{x) is electric field. We simu- 
late transport for 10'* electrons at each value of injector- 
detector voltage drop and the arrival time at the detec- 
tor for each is recorded. The distribution of arrival times 
/(<:, V) is constructed from a histogram of this data and 
are used to calculate the expected output due to spin pre- 
cession in a perpendicular magnetic field (Hanle effect): 

/•OC 

Ic2{V)^ / f{t,V){cos^9cosut + sm^e)e-^'''^idt. 
Jo 

(1) 

where 9 is the angle between the injector/detector mag- 
netization and the device plane, Tsf is effective spin life- 
time and the spin precession angular frequency to = 
g^BB/h. The tilting angle 9 is caused by the external 
magnetic field partially overcoming the finite geometric 
anisotropy of the magnetic thin filmsi^ and is necessary 
to correctly model the experimental results. As in Rcf.^^., 
we use 9 = sm~^{B/15.5kOe). 

In addition, the final spin polarization after transport 
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FIG. 7: Comparison of experimental and Monte-Carlo simu- 
lated voltage dependence of (a) Hanle peak full width at half 
maximum; and (b) spin polarization V = {Ic2 ~ ^02)/ (^C2 + 
Ic2) 152K. The simulation shown in (a) is insensitive to a 
choice of Tsf over the range 10-lOOns. 



can be calculated from 



r{v) 



(2) 



VI. MODEL RESULTS 

As can be seen in Fig. IH an electric field opposing 
transport to the detector is present at low voltage. Elec- 
trons must therefore diffuse against this electric field to 
escape the confining potential in the bulk of the Si trans- 
port layer. Under these conditions of diffusion-dominated 
transport, the arrival-time distribution has a very wide 
exponential shape with average transit time of approx- 
imately 500ns, as shown in Fig. [5] (a). Although the 
width of the distribution can be reduced significantly by 
increasing the voltage drop to the point where the Si 
transport layer is fully depleted as shown in Fig. O (b), 
the confining electric field remains and the exponential 



shape is maintained. This indicates that diffusion is still 
strong. 

For sufficiently high voltage drops, the potential en- 
ergy minimum is annihilated by the detector boundary 
as indicated in Fig. 3] and drift-dominated transport oc- 
curs. This is reflected in the gaussian-like shape of the 
distribution in Fig. [5] (c) for a voltage drop of 6V. At 
this voltage, the average transit time is only approxi- 
mately 50ps, consistent with the analysis of experimental 
Hanle effect measurements. Therefore, as a result of our 
Monte-Carlo modeling we see that the average electron 
transit time in our doped Si spin transport devices can 
be controlled over approximately 4 orders of magnitude 
by changing the injector-detector voltage drop by only 
several volts (from 3V to 6V). 

Using Eq. (jT]), we simulate the Hanle effect in our 
devices using Tgf — 100ns (choice of this value will be 
discussed later). Fig. [6] shows Hanle effect simulations 
for voltages corresponding to the same Vci values as in 
Fig. [3] (again, a shift of 0.3V due to a portion of the 
emitter bias dropping across the resistive tunnel junction 
base^S is accounted for to make a direct comparison) in 
wide agreement to those experiments. In particular, the 
qualitative shape and precession minima positions are 
well modeled. 

The most salient feature of the Hanle effect simulation 
is the magnetic-field width of the central (zero precession- 
angle) peak, plotted as a function of injector-detector 
voltage drop in Fig. [TJ^a) and compared to the experimen- 
tal values. Note that the width is constant for voltages 
greater than 6V (due to drift velocity saturation at high 
electric field in Si), and the presence of a threshold near 
5V (due to appreciable lowering of the confining potential 
barrier at the detector side of the transit layer once full- 
depletion occurs at approximately that voltage). This 
sudden collapse of the Hanle peak width is not seen in 
the voltage dependence of spin precession measurements 
using undoped drift-dominated spin transport devices. 



In Fig. [7] (b), we show the voltage dependence of the 
measured spin polarization V = {Iq2 ~ ^C2 ) / i^C2 + ^C2 ) 
from experimental data using in-plane magnetic field 
spectroscopy as described in Fig. [21 Again, a threshold 
is seen in the experimental data. However, the position 
of the spin polarization threshold in Fig. [7](b) near 3.5V 
is at much smaller bias voltage as compared to the Hanle 
width collapse threshold near 5V shown in Fig[7I^a). This 
indicates that the electrons maintain their spin despite 
a long dwell time which causes strong spin dephasing in 
the confining conduction band potential minimum at low 
voltages. Comparing this behavior to the model results 
from Eq. Q with different values of spin lifetime Tsf 
shows that this discrepancy in threshold position is con- 
sistent with a long spin lifetime of 10-lOOns. This can be 
compared to a spin lifetime of approximately 73ns mea- 
sured in undoped Si at the same temperature using a 
different techniqueJ^ 
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VII. CONCLUSION 

In summary, we have demonstrated spin transport 
through n-typc doped Si. Using a Montc-Carlo algo- 
rithm to model drift and diffusion, we simulated elec- 
tron transport through the inhomogeneous internal elec- 
tric field and make quantitative comparisons to experi- 
mental values of spin polarization and Hanle peak width 
without any free fitting parameters. Analysis of the 
arrival-time distribution indicates that in doped trans- 
port layers, the spin-polarized electron transit time can 



be controlled over several orders of magnitude with ap- 
plied voltage. The resulting non-ohmic behavior seen 
here is in contrast to spin transport measurements using 
undoped silicon transport layers, and is expected to influ- 
ence future semiconductor spintronic device designs uti- 
lizing current-sensing spin detection methods in n-type 
doped semiconductors. 
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